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The structure of bundles of single- walled carbon nanotubes (SWNT) has been refined by Rietveld 
analysis using neutron and X-ray powder diffraction data. Based on previous simulation studies of 
powder diffraction data of SWNT and standard Rietveld analyses, we have developed a pattern fit 
technique for SWNT which provides precise structure parameters. We also show that the present 
technique can be used with the maximum entropy method (MEM), which is complementary to 
the Rietveld analysis. Using the neutron diffraction data of pristine SWNT, we have successfully 
reconstructed the density of carbon nuclei and zero density in the inner cavity of SWNT by MEM. 

PACS numbers: 



I. INTRODUCTION 

In recent years, adsorption of atoms and molecules in 
bundles of single-walled carbon nanotubes^'^ (SWNT) 
is attracting interests for both fundamental and ap- 
plied scientists;^ ranging from investigations of quasi- 
one-dimensional phases to applications of gas storage.^ 
SWNTs are self-organized into a bundle with a two di- 
mensional (2D) triangular lattice structurei^ In this bun- 
dle there are three main adsorption sites: an inner cavity 
of SWNT (tube or t), an interstitial channel (IC) be- 
tween three SWNTs, a groove separating two SWNTs on 
the external surface of the bundle [see Fig. ^a)].^ For 
a small adsorbate, gas uptake may occur progressively 
on the three sites and finally covering the outside of the 
bundle as pressure is increased. Experimentally, gas up- 
takes can be studied by measuring adsorption isotherms. 
A more direct method of studying the adsorption sites, 
particularly useful for the IC and tube sites, is a diffrac- 
tion experiment and its structure analysis of the 2D tri- 
angular lattice. Since the diffraction pattern is not af- 
fected by uptakes by impurities, such as graphite and 
nano-particles of graphite, the diffraction study provides 
complementary information to the adsorption isotherm. 

Techniques of X-ray and neutron diffraction were used 
to demonstrate the triangular lattice of the bundles of 
SWNT,^ and to measure the diameter distribution of 
SWNT.^-^ Powder diffraction patterns by bundles of 
SWNT with several adsorbates have been studied to de- 
termine the adsorption sites for O2 and Nt^^ CDif^'^ 
Ar^ D2rS and I„i^ As pointed out in Ref. I3, intensity 
of the (10) reflection increases or decreases depending on 
adsorption at the IC or tube sites, respectively, being 
used as a simple experimental distinction. 

When internal structures of adsorbed atoms or 



molecules in the tube site are studied, entire powder 
diffraction patterns have to be analyzed more carefully. 
Among a few examples of this kind, Cgo and H2O 
molecules in SWNT were shown to form interesting one- 
dimensional lattices in the inner cavity channel using 
X-ray powder diffraction&i^ti^ Adsorption of iodine in 
SWNT and its clustering in the tube and IC sites were 
studied using both X-ray and neutron powder diffrac- 
tion.^^ However the analyses of these powder diffraction 
patterns are in a primitive stage in comparison with re- 
cent elaborate techniques of the Rietveld analyses for 
three dimensional (3D) crystals»i& The major difficulty 
for SWNT is brought about by randomness existing in 
bundles which is not controlled even by recent prepara- 
tion methods. A bundle consists of different SWNTs with 
numerous chiralities, which are only characterized by a 
radius distribution with an average i? ~ 7 A and a width 
Ai? - 1 A [fuU width at half maximum (FWHM)]. In ad- 
dition, radii of bundles in a sample vary in a typical range 
from ~ 2R to ~ 20R. Small diameter bundles give rise 
to large broadening of Bragg peaks, which complicates 
the pattern fit method. 

The first purpose of the present work is to extend the 
previous simulation-type studies of powder diffraction for 
SWNTiSiSiH to a quantitative Rietveld fit technique. We 
tried to provide a robust basis to this Rietveld analysis, 
which will facilitate to extract maximal information of 
the structure of adsorbates. In discussing this analy- 
sis in detail, it becomes clear that the powder diffrac- 
tion patterns are more informative than previously ex- 
pected. This leads us to another extension of the simula- 
tion studies, a possibility to exploit a Maximum entropy 
mcthodiiiS (MEM) for SWNT. Thus our second pur- 
pose is to formulate the MEM for SWNT, and to show 
its simplest application. This MEM can reconstruct the 
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density of carbon and adsorbates at the tube site without 
assuming any structure models, which can be a comple- 
mentary method to the Rietveld analysis, where a model 
structure is assumed from the beginning. In this paper, 
we restrict ourselves to methodologies and to analyzing 
powder diffraction data of pristine SWNT samples with 
no nominal adsorption, showing simple demonstrations 
of the Rietveld analysis and MEM. 



bundle of SWNTs 



II. 



EXPERIMENTAL 



The raw soot of SWNT was prepared by laser vaporiza- 
tion- of a carbon rod including Ni and Co catalysts, and 
then purified as previously reportediiS The average diam- 
eter of the SWNT used was 13.6 A close to 13.7 A of the 
(10,10) SWNT. To remove adsorbed gases, the sample 
of SWNT was heated at ~ 100°C in an evacuated alu- 
minum container. A neutron diffraction experiment was 
performed using the triple-axis spectrometer 4G-GPTAS 
installed at JRR-3M JAERI (Tokai). The spectrome- 
ter was operated in a triple-axis configuration with coUi- 
mations 30'-80'-80'-open. Neutrons of the wave length 
A = 2.44 A were selected by the pyrolytic-graphite (002) 
monochromator and analyzer. Higher-order neutrons 
were removed by the pyrolytic-graphite filter. 



III. 



STRUCTURE MODEL 



A simple structure mode&Sii of a bundle of SWNTs for 
analyses of diffraction experiments is the 2D triangular 
lattice of cylinders which is shown in Fig.^a). In what 
follows, we use the terminology of neutron diffraction for 
clarity, which can be easily translated into that of X-ray 
diffraction. The cylinder represents the average nuclear 
density of carbon nuclei of SWNT in a bundle. The av- 
erage density is regarded to have cylindrical symmetry, 
where atomic network structures of the graphene sheets 
are ignored, because SWNTs in the bundles are assumed 
to have various chiralities and to be self-organized with- 
out azimuthal correlation owing to the weak coupling by 
the van der Waals force 

We assume that the nuclear density weighted by the 
scattering length of carbon nuclei 6c is well approximated 
by a Gaussian distribution with an average radius R and 
a standard deviation aji (FWHM is Ai? = VSln 2a r), 
which is illustrated in Fig.Q^c). The density around the 
tube site pt{r) — pcir) is given by 



Pcir) = bcac 



2'KaR 



exp 



1 (r-R 

2 [ cjR 



(1) 



where r and trc stand for the radius of the cylindrical 
coordinates and the surface density of carbon nuclei of 
the graphene sheet, respectively. This uniform cylinder 
model (Ai? — > is usually assumed) on the triangular 




ciirxTT interstitial . .\ 
^^^^ channel (IC) , 
\v . cavity (t) 

unit cell 




FIG. 1: (Color online) Structure model of a bundle of SWNTs 
for analyses of powder diffraction experiments. SWNTs are 
approximated by uniform cylinders, in which the network 
structure of the graphene sheet is neglected. They crystal- 
lize into the 2D triangular lattice as shown in (a), where the 
interstitial channel (IC) and the inside cavity of SWNT (t) 
are shown by dashed lines. The unit cell of the triangular lat- 
tice is the region surrounded by the dashed line in (b), where 
the density of the six ICs are multiplied by a weight factor 
of i in calculating the structure factor. The nuclear densities 
weighted by scattering lengths of nuclei at the tube and IC 
sites, Eqs. Q and ||5J with parameters of the Rietveld fit re- 
sult in Sec. II VI are illustrated in (c), where r represents radius 
from the tube or IC centers. 



lattice was proposed to account for images of transmis- 
sion electron microscope and X-ray powder diffraction 
patterns:- The fact that the atomic network can be ne- 
glected was justified by numerical calculations of powder 
diffraction patterns in a low wave-vector range of Q < 2 
A"-'^, which were performed for SWNT with single chi- 
rality.^ 

The nano-scale channels in the bundle denoted by the 
tube (t) and IC sites are illustrated in Figs. H^a) and 
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(b). For calculations of diffraction patterns, adsorbed 
atoms in the tube site can be treated by an additional 
contribution to pt{r) — pci'i') + Pads('')- These atoms 
can show various structures and states in the wide inside 
space of SWNT, which should be studied by diffraction 
experimentsiSii^^ On the other hand, adsorbates in the 
IC site are localized in the narrow space surrounded by 
the three SWNTs. For the analysis of powder diffraction 
in the range Q < 2 A~^, we can safely assume that the 
nuclear density is well approximated by a Gaussian func- 
tion centered at the IC site with a small width an^^^ (half 
width at half maximum is Ai?ic — V2 In 2a i^^^ = 1 A). 
This nuclear density weighted by the scattering length of 
nuclei in the IC site bic, which is shown in Fig. ^c), is 
written by 



pic(r) = feicpicexp 



1 / r 
'2 \^ 



(2) 



where r and pic are the radius from the IC center and 
the peak nuclear density, respectively. 

The unit cell of the 2D triangular lattice, which is 
used in calculating the structure factor, is depicted by 
the dashed line in Fig. ^b). It should be noted that, to 
keep the six-fold symmetry of the structure factor as a 
function of wave vector, we choose a non-standard unit 
cell containing one tube site and six IC sites, the latter 
of which is weighted b y ^ f or the structure factor calcu- 
lation (cf. Eqs. IIA.6|) . (|A.7|) . and IjA.SI) in the appendix). 

The average radius R of SWNT in the bundle de- 
pends on its condition of the preparation, e.g. the lo- 
cal temperature of the furnace, implying that R takes 
a number of values for bundles in a sample. Following 
Refs.HandH we assume that the probability density of 
R is approximated by a Gaussian distribution with an 
average i?av and a standard deviation a^^^ (FWHM is 



VS In 2(Tfl^^ ) , expressed by 



piR) 



1 



27rCT 



■ exp 



1 f R — R„ 



(3) 



We note that the powder diffraction intensity by a bun- 
dle, which is calculated using a single value of i?, has to 
be further averaged over the probability density. In this 
average, R dependence of the intertube gap g — a — 2R 
[see Fig.njb)], where a is the lattice constant of the tri- 
angular lattice, can be neglectedi^ 



IV. RIETVELD ANALYSIS 

In standard Rietveld analyses!^ of powder diffraction 
data of 3D crystals, the intensity of the diffracted beam 
at a scattering angle 29 is fit to 



I{29) ^sA{20) mK\FK\^L{29K)H'20 - 29k] 



K 



(4) 



yb(20), 



where rriK, Fk and L{29k) = l/(sin6'if sin 26*/^ ) are the 
multiplicity, the structure factor, and the Lorentz factor, 
respectively, of the K-th. reflection with the scattering an- 
gle 29k] s, A{29) and yh{29) represent the scale factor, 
the absorption factor, and the background, respectively. 
The profile function ${26* — 29k), which is normalized by 
/ $(26* — 29K)d{29) = 1, shows if-dependent broadening 
due to the instrumental resolution and the sample crys- 
tallinity. The functional form of is often approximated 
by the pseudo-Voigt function, which is further modified 
to an asymmetric form^SiSS 

For a bundle of SWNTs, the scattering cross-section 
due to the structure model in the previous section, which 
is derived in the appendix [cf. Eq. (|A.ll|l ] , is given by 



da 



iQ) 



V3a^ ^ 



(5) 



y-\nQ±)\\ 



where and Nt stand for the length of SWNT [see 
Fig. n^a)] and the number of SWNTs in the bundle, 
respectively. The wave vector Q = QzG-z + Q± is de- 
composed into the z component Qz parallel to the di- 
rection along SWNT, and its perpendicular components 
Qj^ = ha* + kb*. The structure factor F{Q±) is defined 
by Eq. HA.4() . and the function S {Q± — G) is the broad- 
ened 2D 5 function peaked at a reciprocal lattice point G 
of the triangular lattice. Simulations of powder diffrac- 
tion by bundles of SWNTs have been performed to study 
X-ray and neutron diffraction datar'^i6i7i8iiQiiiii2,i3,i4,i5,2i 

As noticed by Thess et al.,^ an important difference of 
Eq. (jSJ from that of the 3D crystals is the broadening of 
the 2D 5 function by the finite size effect of the bundle ra- 
dius. This broadening of the Bragg peaks is much larger 
than instrumental resolutions of powder diffractometers. 

Powder averaged intensity of Eq. (0 , corresponding to 
the first term of Eq. Q, is given by [cf. Eq. HA.14|) ] 



i{Q) c^sA{Q)Y,^k\F{QGk)?L{Q) 



K 



(6) 



x<t>{Q-\GK\), 



where F{QGk) is an approximation of the powder av- 
eraged structure factor [cf. Eqs. (jA.ll|l . (jA.12|l . and 
(jA.13|) ]. and (f){Q — |C?a'|) is the normalized profile 
function, which is further modified to an asymmetric 
form>iS*2S The use the wave vector Q = \Q\ instead of 
scattering angle 29 is the previous convention, ^■^■'^ where 
the Lorentz factor becomes L{Q) — l/sin'^6' or oc 
By further taking the average of Eq. ^ over the distri- 
bution of the radius R of SWNT, Eq. (|5J), we arrive at 
an expression 



/(Q) = (/(Q)>i^ + ^/b(Q) 



(HQ)) 



R 



piR)i{Q)dR, 



(7) 



which is used as the calculated intensity for the present 
Rietveld fit. 
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Neutron powder diffraction pattern of pristine SWNT 
with no nominal adsorbates measured at 300 K is shown 
in Fig. 12 The instrumental resolution with the present 
condition is not negligibly small compared to the broad- 
ening by the finite size effect. This additional resolu- 
tion correction was performed by convoluting {I(Q)) r of 
Eq. O with the Gaussian resolution function:-?^, The pro- 
file function (/){Q — |Gi<-|) of Eq. © was determined by 
the similar Rietveld fit based on Eq. using the pre- 
vious X-ray diffraction datck^ of pristine SWNT. These 
data were measured using a synchrotron X-ray powder 
diffractometer with far higher Q resolution, where no res- 
olution convolution was made. The profile function could 
be well approximated by superposition of four Gaussian 
functions. The resulting Rietveld analysis is plotted in 
the inset of Fig. [3 Using this profile function, the Ri- 
etveld analysis of the neutron powder diffraction was per- 
formed. The background yh{Q) of Eq. iQ was chosen to 
be a polynomial in Q"^, X]n=o ''"'5^"5 with additional 
three Gaussian functions centered around Q ~ 1.8 A^^, 
where impurity peaks by graphite and nano-particles of 
graphite exist. The fitted curve shown in Fig.|2]is in good 
agreement with the observed data. 

The fitted parameters of the model structure are: 
i?av = 6.78 ± 0.02 A, Ai?av = 0.56 ± 0.04 A (FWHM), 
Ai? = 1.32 ± 0.08 A (FWHM), g = 3.34 ± 0.04 A, 
Ai?ic = 1 A (fixed), and hcpic/ ibcoc) = 0.21 ± 0.08 
A~^, for the neutron diffraction data. Those for the 
X-ray diffraction data are: 

i?av = 6.715 ± 0.005 A, 
Ai?av = 0.68 ± 0.01 A (FWHM), AR = 1.26 ± 0.03 
A (FWHM), g = 3.23 ± 0.01 A, Ai?ic = 1 A (fixed), 
and Po,ic/o'o,c = 0.11 ± 0.04 A~^, where Pcjc and CTc,c 
are the peak electron density at the IC site and the sur- 
face electron density of the graphene sheet, respectively. 
These parameters of the X-ray data are consistent with 
the previous work.^ They also agree reasonably well with 
the parameters of the neutron data, considering the slight 
sample difference. We conclude that the Rietveld fittings 
have been successfully performed for the neutron and X- 
ray powder diffraction data. 



neutron powder diffraction JRR3M-4G ki=2.57 30'-80'-80'-open 

7000 



V. 



MAXIMUM ENTROPY METHOD 



The structural information of the unit cell is contained 
in the structure factor F{Q). In a Rietveld analysis of an 
ordinary 3D crystal, structure factors are observed at a 
number of discrete wave vectors G. On the other hand, 
there are only 10 refiections in the present neutron pow- 
der diffraction data of SWNT (see Fig. [2J|. From this 
small number, one may think that it is very difficult to 
determine many structural parameters. However, this is 
not simply the case, owing to the other aspect of the 
important difference, i.e., the large width of the profile 
function, showing substantial overlaps with neighboring 
reflections, which is easily seen in the range 0.4 < Q < 0.8 
A-iofFig.|21 This overlap suggests a possibility to make 
use of the structure factor in a continuous Q range. If the 
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FIG. 2; (Color online) Neutron diffraction pattern of SWNTs 
at 300 K. Observed and calculated intensities are plotted by 
closed circle and dotted lines, respectively. Background is 
plotted by dashed line. Vertical bars stand for Bragg peak 
positions with a — 2i?av + <?. X-ray powder diffraction of 
the previous work*" and its Rietveld analysis are shown in the 
inset. 

observer range of 0.2 < Q < 2.1 A-i is fully available, 
one may think that the powder diffraction data contain 
more information than the simple expectation. As no- 
ticed in Ref . '21', this is the reason why the parameter i?av 
is precisely determined by powder diffraction of SWNT. 
In the following, we discuss a method which takes advan- 
tage of the availability of the continuous Q range. 

Adsorbates in the IC site, which are localized in the 
small space, cannot have internal structures except for 
those in the direction along SWNT. Thus we assume that 
the IC part of the structure factor has only one parameter 
Pic, which can be easily determined. Then remaining 
problem is the structure at the tube site, SWNT and its 
inside, which is expressed by Ptif)- The density ptir) is 
related to the tube part of the structure factor by [cf. 
Eq. iXflll ] 



F,{Q) = 2tt rJo{Qr)ptir)dr, 



(8) 



where Jo{Qr) is the cylindrical Bessel function. Thus 
the problem is reduced to how to determine pt(r) from 
Ft{Q). Since Eq. (O is the Fourier-Bessel transform, a 
solution is to use the inverse transform 



1 f^- 



QJo{Qr)Ft{Q)dQ. 



(9) 



However an application of this inversion to the neutron 
diffraction data failed, because the inverted ptir) showed 
a spurious oscillation, which is a cut-off effect due to in- 
sufficient Q range. A well known method to avoid this 
cut-off effect of Fourier transform is MEMjiii^i which 
is being successfully used in single-crystal and powder 
diffraction analyses 
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In general, MEM is based on a linear relation between 
observed and calculated quantities (data and image in the 
MEM terminology^). This relation is Eq. © for Ft{Q) 
(data) and pt{r) (image). To apply a general purpose 
MEM algorithmj^^ a discrete version of the linear relation 

FtiQn) = 27r^ Jo{Qnrm)rm^rmPtirm) (10) 

m 

was used as a basis for the present MEM. The first step 
of the MEM is to determine observed values of Ft{Qn)- 
We calculated these by a Rietveld fit, where Ft{Qn) are 
fitting parameters and intermediate Ft{Q) is evaluated 
by a linear interpolation. In this fitting the other param- 
eters: the background parameters of yh{Q), s, &icPic> 
Ai?ic, 5, Rav, and Ai?av, were fixed to the values of 
the previous section. It should be noted that the three 
parameters, g, i?av, and A_Rav, and the SWNT radius 
R still have meaning in connection with the lattice con- 
stant a — 2R + g. Thus the fitting parameters Ft{Qn) 
depend on R, and we refer to this R dependence explic- 
itly by [Ft{Qn)]ji- In performing the average over R [cf. 
Eq. ITJl], we assumed R dependence 

{FtiQn)]n = ^F^'^ (^Q?j, (11) 

which is a good approximation for empty SWNT. The 
Rietveld fit was carried out and resulting parameters 
[^t(Qn)]fl=/s are shown in Fig. |2Ia) together with that 
of the Rietveld analysis in the previous section. 

The second step of the present MEM is to reconstruct 
Ptirm) from [Ft{Qn)\ii=R^^- carried out this using the 
general purpose MEM algorithm.^'' The reconstructed 
Pt{i'm) is plotted in Fig.|3Jb), where the Gaussian func- 
tion Eq. with R ~ i?av determined by the Rietveld 
analysis in the previous section is also shown. One can 
see that the Rietveld and MEM results agree surprisingly 
well. The structure factors [Ft{Qn)]]:j-ji were calcu- 
lated using the MEM reconstructed pt{rm) and Eq. (|10|l . 
and are plotted in Fig. |2|;a). They agree well with the 
observed [Ft{Qn)]fi^fi apart from small discrepancy 
around Q ~ 0.55 A~^, implying that the present MEM 
gives self-consistent results. Therefore we conclude that 
the continuous Q data of structure factor have been ex- 
cellently analyzed by the present MEM formulation. 



VI. DISCUSSION 

The advantage of MEMiiS^ is providing a way to 
choose a maximally non-committal nuclear density with- 
out any structure models among a series of nuclear den- 
sities consistent with experimental data. Thus MEM 
is a complementary method to the Rietveld analysis, in 
which an empirical structure model is assumed from the 
beginning. Since the formulation of MEM discussed in 




FIG. 3: (Color online) (a) Structure factors of the tube site 
Ft(Q) with R = _Rav are plotted as a function of wave vector. 
The full line is the result of the Rietveld analysis in Sec. IIVI 
Closed circle represents the observed [Ft{Qn)]n^j^ deter- 
mined by the fit described in the text. Open square is the 
calculated [-Ft(Q„)]^^^^^ using Eq. HUH and pt{rm) obtained 
by MEM. (b) Nuclear densities at the tube site obtained by 
the Rietveld analysis in Sec. |T3 with R = i?av and MEM. 



Sec. partly relies on the result of the Rietveld analy- 
sis of Sec. IIVI we have to check consistency between the 
MEM reconstructed nuclear density and that of the em- 
pirical structure model. Fig.|3^b) definitely warrants this 
consistency for the pristine SWNT sample. For general 
cases, where atoms are adsorbed inside SWNT, the con- 
sistency is not necessarily obtained after the first trial 
of MEM. If the first consistency check is negative, we 
have to adjust the structure model considering the first 
MEM reconstructed nuclear density and perform second 
Rietveld and MEM analyses. Depending on the second 
consistency check, we may arrive at a final answer or 
continue this iterative procedure. This kind of iterative 
use of the Rietveld and MEM analyses have been skill- 
fully used for powder diffraction of 3D crystals>iSi24 Ap- 
plications of the present Rietveld and MEM analyses to 
SWNT with adsorbates in the tube site are awaited. 
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VII. CONCLUSIONS 

We have developed a pattern fit technique, Rietveld 
analysis, for powder diffraction data of bundles of 
SWNTs, by extending previous simulation studies. Neu- 
tron and X-ray powder diffraction data of pristine SWNT 
samples have been successfully analyzed by this Rietveld 
analysis. The high fitting quality indicates that the anal- 
ysis will be useful to obtain structure parameters of ad- 
sorbates in SWNT samples. We also have exploited a 
MEM for reconstructing nuclear or electron density at 
the tube site with the use of the pattern fit technique. We 
have applied this MEM to the neutron diffraction data 
and reconstructed the nuclear density of SWNT consis- 
tent with the Rietveld result. 



where the summation is taken over 2D triangular lattice 
positions. The structure factor F{Qj_) is given by 



F{Q±)^ ( exp(iQj 

J unit cell 



(A.4) 



where the unit cell is shown in Fig. ^b). We note that 
although the 5 function in Eq. IjA.Sp is sharp along the 
Qz direction, the lattice sum peaks at reciprocal lattice 
points G with a finite width in the Q± space. ^ The lat- 
ter would become 2D 5 function only in the limit of the 
infinitely large radius of the bundle, as 



^exp(zQ_L • ■r_L,t) 



^Y.5{Q^-G), (A.5) 
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APPENDIX: SCATTERING CROSS-SECTION 



where Nt{-^ oo) is the number of SWNTs. 

The finite width of the Bragg-like scattering in 
Eq. ljA.3(l by small A't requires us to carefully treat the 
structure factor of Eq. (|A.4p . because F{Q±) should 
be evaluated for continuously varying Q±. In view of 
the rope shape of TEM images of bundles,^ which has 
quasi-six-fold symmetry, it is probably adequate to keep 
the six- fold symmetry of F{Q±) for continuously vary- 
ing Q±. It can be retained for the unit cell shown in 
Fig. n^b) , explicitly using the nuclear density 



The neutron scattering cross-section by a bundle of 
SWNTs is calculated using the general formulck^^ 



da 



(Q) 



J exp{iQ ■ r)p{r)dr 



(A.l) 



where p(r) stands for the nuclear density weighted by the 
scattering length of each nucleus. On the basis of this 
equation one can discuss the Bragg scattering, scattering 
from a single SWNT, cluster of atoms or molecules, small 
angle scattering, etc. 

For the Bragg scattering due to the model structure 
shown in Fig. ^ the nuclear density is expressed by 



p{r ^r^ + ze^) = ^ Pt{\r± - r_L,t|) 

+ X! P^c{\r± - r±jc\)- 



(A.2) 



The first term is the nuclear density of SWNT and adsor- 
bates at the tube site r_L,t = naa + ni,b. The second term 
is the nuclear density of adsorbates at the IC site rj__ic — 
r_L,t±(ia-fJb) [or± (|a+ ifo) ,± (ia- ib)]. By in- 
serting Eq. ljA.2|) in Eq. IjA.ip and performing the inte- 
gration along the z axis in < z < with large Lz > 1 
p,m,i one has 



^ (Q) =27rLz6{Qz 



X \F{Qj 



^ exp(iQ_L ■ r_L,t) 



(A.3) 



1 ® 

p{r±) = Pt{\r±\) + 1^"^ pic{\r± - r±^ic,i 



(A.6) 



in Eq. IIA.4|) . where rj^jc,,; is the surrounding six IC sites 
around the tube site r±^t = 0. Substituting Eq. l|A.6p for 
Eq. HA.4|) . the structure factor becomes 



F{Q±) = Ft{\Q^\) + 9{Q±)Fic{\Qj 



(A.7) 



where 
g{Q± = ha* 



kb*) = -< cos 
^ 3 



^{h + 2k) 



27r 

— {2h + k) 



(A, 



The two terms in Eq. HA.7|I are the contribution from the 
tube site 



Ft{\Q±\) = J cxp{iQ^-r^)pt{\r^\)dr^ 



(A.9) 



2tt / rJo{Q±r)pt{r)dr, 

JQ 



and that from the IC site 



Fic{\Q±\) = y"cxp(iQ_L ■ rj_)pic{\rj_\)drj_ 



(A.IO) 



27r / rJo{Q±r)picir)dr. 
Jo 
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Considering Eqs. ljA.3|l and ljA.5|) . the neutron scatter- 
ing cross-section due to the Bragg-like scattering of the 
2D triangular lattice is summarized as 



da 



x\nQ±)\\ 



6{Q,)S{Q^-G) 



(A.ll) 



where the function 6 {Q± — G) is a broadened 2D 6 func- 
tion normalized by J 5 {Q± — G) dQ± = 1. 

To calculate intensity of powder diffraction, this cross- 
section is further averaged over all orientation of the bun- 
dle, or equivalently over all direction of Qr^ i.e., 



da 

dn 

1 



(Q) 



Ql^Q, powder av. 



47rg2 



da 



(A.12) 



iQ)dSQ, 



where dSQ is the infinitesimal area on the sphere of ra- 
dius Q. Since the numerical evaluation of this powder 
average took too long time to perform the Rietveld fit- 



ting, we made an ad hoc approximation given by 

1 



\G\ 



G 



(A.13) 



V3a2 
F{QG)\ 



where a profile function (f>{Q — G) normalized by J ( 
G)dQ = f is used, and G — G/G. We checked numer- 
ically that Eq. (|A.13|I is a good approximation for the 
fits of the present experimental data. By using the six- 
fold symmetry of the structure factor, wc finally write 
Eq. IjA.f 3(1 in a form appropriate for the powder diffrac- 
tion 



i{Q) ^L,N,^ Y.'^Km~\GK\)^ 



(A.14) 



X F{QGk)\ , 



where mn is the multiplicity. 
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